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Abstract. The closed form integral representation for the projection onto the 
subspace spanned by bound states of the two-body Coulomb Hamiltonian is obtained. 
The projection operator onto the n 2 dimensional subspace corresponding to the rt-th 
eigenvalue in the Coulomb discrete spectrum is also represented as the combination of 
Laguerre polynomials of n-th and (n — l)-th order. The latter allows us to derive an 
analog of the Christoffel-Darboux summation formula for the Laguerre polynomials. 
The representations obtained are believed to be helpful in solving the breakup problem 
in a system of three charged particles where the correct treatment of infinitely many 
bound states in two body subsystems is one of the most difficult technical problems. 
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1. Introduction 

The two-body Coulomb problem is perhaps the most famous problem of quantum 
mechanics formulated on the basis of the Schrodinger wave equation pQ. First solved 
analytically for the bound states pQ, it was then solved for scattering states j2j, as 
well as for the Green's function [Hj, [I], (Sj- The situation is quite different for the 
systems involving three and more charged particles. There is no analytic solution to the 
Schrodinger equation in this case and the problem exhibits a great complexity especially 
if the ionization process is energetically permitted. In this case, the infinitely many open 
excitation channels are lead to extremely complicated behavior of the wave function, 
which asymptotically possesses infinitely many terms. Although definite progress has 
been made in the practical numerical solution of the three charged particles problem 
above the disintegration threshold fH| by methods avoiding the explicit use of the wave 
function asymptotics, the theoretical status of the few-body Coulomb problem is still 
unsatisfactory in many respects. Among the very extensive literature devoted to the 
few-body Coulomb problem, the following works [7j, [H] and refer to the theoretical 
aspects of the problem. 

Recently we have presented a new method of handling the Coulomb potentials in the 
few-body Hamiltonian with the help of the Coulomb- Fourier transform (CFT) JU]. The 
method allows us to exclude the long-range Coulomb interaction from the Hamiltonian 
by a specially constructed unitary transformation. This method was proven useful for 
repulsive Coulomb interactions. In the case of attraction, the analytic closed form 
representation for the projections onto the Coulomb bound-states subspace may lead to 
a substantial simplification of the CFT method machinery. Although the latter was the 
primary goal, it was found that the representations obtained, being quite general, are 
not well addressed in the literature, and this has stimulated this publication. 

The paper is organized from three sections. After the introduction, in the section 
two we derive representations for projection operators and consider some particular 
cases. The third section concludes the paper. Throughout the paper the bold letters, 
e.g. r, r', are used for vectors and not bold for their magnitudes, e.g. r = |r|. The unit 
vector associated with r will be denoted by r = r/r. 

2. Representations for projections 

We consider the infinite dimensional projection operator 



The operators V n are the orthogonal projections onto the n 2 -dimensional subspaces 
spanned by the two-body Coulomb bound-states (r\ip n [ m ) = ip n i m (r). V n have kernels 
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where the normalized Coulomb bound-state wave functions are chosen in the form 
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Vw(r) = a 3/2 — , 



n- 



f^^Fj^^™ (3, 

Here a = ^ r Z 1 Z 2 >0, functions L$ and Y™ are the generalized Laguerre polynomials 
and spherical harmonics as they are defined in (TTj, respectively. The wave-function 
■^nZmfr) obeys the Schrodinger equation 

2a 

(H - E n )^ nml {r) = (-A r E n )iP nml (r) = (4) 

with E n = —a 2 /n 2 and n positive integer. 

In order to work out the representations for Vd, n we are seeking, let us begin with 
the standard formula for the projection V n as the residue of the Green's function 

p n (r, r') = / £ c (r,r',CR (5) 

Z7T1 JC En 

with the contour Cs n encircling the point E n in positive direction in the ( complex 
plane. Then, using the Hostler [3] representation for Q c (r, r', Q 

where u = a/y/( and s± = r + r' ± |r — r'| and evaluating the residue, we arrive at the 
following expression for V n (r, r') 

V n (v,v') = -- . (7) 

n 4 ir{s + — s_J 

fas-\ (i) fas + \ (i) fas-\ (as+ 



x 

n J " " V n / V n J \ n 

Alternatively, we can transform (J7|) into the form 

a 2 e-^ s + +s -^ 

VnM = ~ r (8) 

n z ir{s + — s_J 



ray V n / \ n / \ n 

The formulae and (jHJ are the starting points for representations of this paper. 

S.i. Integral representation for projections onto the discrete spectrum subspace 

To proceed with formulae (J7J) and (JHJ), we use the well known expression for Laguerre 
polynomials ^2] i n terms of Bessel functions 

L f_ x (£j = ?L_ n f>+ly-P/2 e v/n J™ dx X P/2~1 ( xe -)" Jf) (2 ^) (9) 

and the following integral representations (3.382.7) ^3] 
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J™dye- nv I (y/2ya(a + bf) 



\ e a{a+b)/2n 
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where Jp{z) and Iq{z) are the Bessel function and the Bessel function of imaginary 
argument, respectively. Introducing these representations into (JJJ) we arrive at the 
following integral for V n (r, r') 

T ^ = Wr^r- (10) 



/oo roo poo / i \ roo f 
dyi dy 2 dx 3 I [J2x 3 a(s + + S-)) / 
-oo J-oo JO V v / JO X\ JO 



00 dxi f°° dx 2 
x 2 



n / ■ \ n 



x [^as + x 2 J (2^as_Xi) Ji{2yjas + x 2 ) — ^ as _x 2 Ji{2^f as _x 2 ) Jq{2^/ as + Xi)] . 
By introducing a new five-dimensional variable X = {x±, x 2 , £3, yx, y 2 } and quantities 



Q = Xl e- Xl+1 x 2 e- X2+1 e- X3 



g-ij/i e~ ly2 
1 - iyx 1 - iy 2 ' 



B(X, a, s + , s_) = I Q \ J2x 3 a(s + + s_ 

x [y/as + x 2 Jo(2y/as-Xi)Ji(2y/as + x 2 ) — y / as~x 2 Ji(2 y /as-X 2 )Jo(2y'as + xi)] 
we rewrite the latter formula in the compact form 

v ^ = jk^-l^f rB(x ' a ' s+ ' s - ) - 

The integration domain Q is defined as 

Q = {X : < Xi < 00, % — 1, 2, 3, -00 <y k < 00, k — 1, 2}. (12) 

Let us notice that the quantity |Q| is bounded on f2 with the only maximum (such that 
|Q| = 1) at the point X = {1,1,0,0,0} , hence everywhere except X the inequality 
|Q| < 1 holds true. 

Before computing the infinite sum (JTJ let us consider the operator 

K: = E v». (13) 

n=Ni 

The kernel of the operator V N \ can easily be computed by using the formula (fTTj) for Pn 
and evaluating of the sum of the geometric progression of Q n terms under the integral 
which yield 

V N N IM = P$(r,r',n) = (14) 
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dX Q JVl (l-Q JV2 - iVl + 1 ^ 



f dX W 1 -g- —J 



Now we are ready to evaluate the limit of P^ 2 (r, r',f2) as N 2 — > oo keeping iVj 
finite. Let us notice that the integral (|14jl converges uniformly at any N 2 and hence for 
any positive e we can find S > such that 

\P$(v,v'MXa,5))\<e (15) 

where Q(Xq,5) is a neighborhood of the point Xq in which the quantity Q reaches its 
maximum, i.e. 

n(X ,6) = {X G Q : \X-X \ < 5}. (16) 

On the rest of the integration domain Q(Xq,5) = fi\0(Xo,5) the inequality |Q| < 1 
holds true and we can take the limit 



P Nl (r,r',Q(X ,S)) = lim P^ 2 (r, r', Q(X Q , 5)) = (17) 

iV2— >oo 



-d(a, a, s+, s_J. (1; 



(27r) 3 s + — s__ 7 xix 2 1 — Q 

n(X ,<5) 

For the integral ()18j) the limit 5 — > is permitted, so that due to the arbitrariness of e 
we get 

V Nl (r,r') = }\mP Nl (r,r'MX ,6)) = (19) 

o— >0 

1 2a 2 f dX Q r 



f dX 



(27r) 3 S + — S_ 7 X\%2 1 — Q 

Now by setting N\ = 1 we arrive at the final result for the projection Vd 



^< r - r » = (2^— / ^ A B(X ■ a - s+ - *->■ (20) 

The formulae (fTTIJ) and ()2())1 are the main results of this subsection. 
2.2. Some particular cases 

In this subsection we consider a particular case of the integral representation (JBJ) which 
leads to an analog of the Christoffel-Darboux formula (22.12.1) jTTj applied to Laguerre 
polynomials. We also evaluate the asymptotics of the projections kernel V n (r, r') in the 
special case when n ^> a(r + r') with the help of representation (JZj). 

It is worthwhile to notice that the formulae ()7I8|) in the particular case when n — 1, 2 

3 

Vl (r,r') = —e- a(r+r,) , (21) 

7T 

p 2 (r, r') = -jie-T (r+p ' ) [4 - 2«(r + r') + a 2 rr(l + f • r')l (22) 

327T L - 1 
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give essentially the same results which can be computed directly from the conventional 
representation ()2I3|) . Let us now consider the case of arbitrary n. From (I2I3|) and (JHJ) 

we get 



1 n ^ (2ar2ar'\ l (n - I - l)\ , (a+1) (2ar\ (2/+1) ( 2ar' s 



(23) 



, E ) ^jf^-i {^) 4SS ^ j (2/ + 1) W • r )• 

Taking a particular case r • f' = 1, setting (2ar)/n = x, (2ar')/n = y and changing the 
summation variable in such a way I = n — m — 1 we get 

(v) 



ar - y 

n-l 



(24) 



1 E {xy)— x m ^ {n m \,} ) L^ n - m) - x \x)L^-^{y). 

n „ 2n — m — 1 ! 



(2r 

m=0 v 

This formula if compared to the Christoffel-Darboux summation formula for Laguerre 
polynomials (22.12.1) [TJJ 

L - l(x)L " (y) ' L " (x)L "" l(y) = i E (25) 

yields the following interesting identity 

n-l 

X! L m {x)L m {y) = (26) 

m=0 

As the last special case we consider the behavior of the projection V n for large value 
of the principal quantum number n. Let us use again the representation Q which we 
rewrite in the form n 

!<2>Mn) = V + V 1/2 ^ /n ^(n,</) 



where F@(n,y) stands for the integral 



FP(n,y)= / dxx^^e-^-^J^^y) 



and introduce it into the formula (jZJ) for Laguerre polynomials. If n ^> a(r + r') then 
the only critical factor under the respective integrals F@(n, as±) is e - n \ x ~ lo & x ) with the 
only critical point x = 1. Evaluating the the integrals F l3 (n,as±) as n — > oo by the 
Laplace method we get for the projections V n (r,r') the following asymptotics 

V n (r,r) oc — — -[ v /s7J (2 v /asT)Ji(2 v /as7) - ^/sIJ 1 (2y/asZ)J (2 y /as + )\. 

Tin ( S-l — S— i 
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3. Conclusion 

The closed form representations are obtained for projections onto the n 2 -dimensional 
subspace spanned by bound-state eigenfunctions of the Coulomb Hamiltonian 
corresponding to the principal quantum number n as well as for the projection onto the 
subspace spanned by all Coulomb bound-states. These representations can be useful for 
solving the few body scattering problem in a system of charged particles for energies 
above the three body disintegration thresholds. The asymptotics computed above for 
the projections V n as n — > oo may lead to drastic simplifications in calculating different 
Coulomb matrix elements between states which are spatially well confined. The analog 
of the Christoffel-Darboux summation formula for Laguerre polynomials which is derived 
as a particular case of the representations for the Coulomb projections can be useful for 
the theory of classical orthogonal polynomials. 
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